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Static Exact Solutions of a Spin Model Exhibiting Glassy Dynamics 

Chiaki Yamaguchi 
Kosugichou 1-359, Kawasaki 211-0063, Japan 

A spin model which exhibits glassy dynamics has been proposed by Newman and Moore. 
This model possesses no randomness for exchange interactions. We propose partial trace 
methods for obtaining static exact solutions of this model under free boundary conditions, 
and show some static exact solutions. We obtain that, from the exact solutions of the specific 
■ heat and the correlation functions, the transition temperature is zero, although the thermal 

average of a spin is zero when the temperature is zero. We investigate the ground states. 
In the ground states, a part of the present results especially disagrees with the part of the 
previous results. We discuss the disagreements. 

§1. Introduction 

Spin glass mo dels ,'3''^ that are known as popular models for glassy systems, have 
randomness for exchange interactions. The mathematical study is hard. In addition, 
there is a problem whether the randomness for exchange interactions is necessity or 
not for glassy properties. Therefore, searching and investigating models, that have 
the properties of glassy systems and have no randomness for exchange interactions, 
can be one of promising methods for understandings of glassy systems. For achieving 
this purpose, a spin model which has no randomness for exchange interactions and 
exhibits glassy dynamics has been proposed by Newman and Moore.'^ For dynam- 
ical features, this model has aging behavior,^ and it falls out of equilibrium at a 
temperature which decreases logarithmically as a function of the cooling time,-' for 
example. Dynamical features of this model are investigated in Refs. [11,13]) ,[5l),EI). 

For static solutions, the solutions have been exactly obtained on lattices of 
length a power of two along at least one dimension under periodic boundary con- 
ditions . '3'® ■'^''^ In other words, this model has a special property. Therefore, this 
model has been solved in a special case that, under periodic boundary conditions, 
the lattice length is a power of two along at least one dimension. There have been 
no solutions on lattices of any length under periodic boundary conditions. In ad- 
dition, there have been no solutions on lattices of any length under free boundary 
conditions. Since the previous method for solving this model uses a special property 
of this model, there is a possibility that a part of the results is different from that on 
lattices of any length, and this issue has not been shown exactly. In this article, by 
solving exactly this model on lattices of any length under free boundary conditions, 
we exactly show that a part of the results is different from the part of the previous 
results. 

We propose partial trace methods for obtaining static exact solutions of this 
model under free boundary conditions, and find some static exact solutions. We 
find the exact solutions under free boundary conditions for the energy, the specific 
heat, the free energy, the entropy, the number of ground-state degeneracy, the ther- 
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Fig. 1. The present model. In this figure, the hnear size of this system, L, is 5, the number of 
spins, A'^, is 24, and the number of plackets, Np, is 16. The plackets are drawn as meshes. 



mal average of a spin, the two-point correlation function, the three-point correlation 
function and the four-point correlation function. We examine the transition temper- 
ature, and investigate the ground states. We discuss the disagreements between a 
part of the present results and the part of the previous results. 

This article is organized as follows. First in ^ the model is explained. In ^ we 
describe the present partial trace methods for solving this model and obtain results. 
In ^ the disagreements between a part of the present results and the part of the 
previous results are discussed. This article is summarized in ^ 

§2. Model 

The Hamiltonian of this model under free boundary conditions, Ti, is given by 

L-l L-1 

'H = j''^''^ax,yCrx+i,yCrx,y+i, (2-1) 

^=1 y=l 

where ax.y is an Ising spin, ax^y = ±1, and x and y are coordinates of the spin. We 
set J > 0. This system is on downward pointing-triangles of the triangular lattice. L 
is the linear size of this system. The number of spins, N, is — 1, and the number 
of plackets, Np, is (L - 1)^. 

Figure [1] shows the present modehl^ In this figure, the linear size of this system, 
L, is 5, the number of spins, A^, is 24, and the number of plackets, Np, is 16. The 



Static Exact Solutions of a Spin Model Exhibiting Glassy Dynamics 



3 



plackets are drawn as meshes. 

Note that this model is distinct from the model of Baxter and Wu,® which has 
interactions on upward-pointing triangles also. Its spin behavior of the model of 
Baxter and Wu is entirely different from this model.'^ 

In the ground state, there is a relation among spins given bjff' 

0'x,y<^x+l,yO'x,y+l = ■ (2"2) 

From this relation, if the boundary spins of system are given, all the spins in the 
ground state are uniquely obtained. 

This model displays glassy behavior under single-spin-flip dynamics. It has been 
pointed out that, in Ref. [1]), this model has aging behavior, and it falls out of 
equilibrium at a temperature which decreases logarithmically as a function of the 
cooling time. See Refs. [T]),[3]),[5]),[6]) for dynamical features of this model. 

The transformed models by a transformation J — t- it J for each three-spin inter- 
action have identical behavior with this model. ^' 

This model under a magnetic field has been also proposed concerned with ther- 
modynamic transitions associated with irregularly ordered ground states.!^ 

§3. Results 

We describe the present partial trace methods for solving this model and obtain 
results. 

We adopt the free boundary conditions as the boundary conditions. Therefore, 
there is a spin which is only on a placket. We perform the integration of the spin at 
first. We define the spin as crj, and define a function A{aj,ak,{cri}), where {ai} is 
the set of all spins except for cjj, aj and cr^. We introduce a variable K and define 
K = (3 J, where /? is the inverse temperature, and /3 = 1/kBT. A:^ is the Boltzmann 
constant, and T is the temperature. Then, we perform the integration of the spin cjj 
as 

ai CFj 

= (e-^ + e"") E E ^(^^•' • (3-1) 

In Eq. (jS-ip . we are able to treat aj and as the spins on the free boundaries. 
Therefore, we are also able to perform the integrations of aj and cjfc by using similar 
integrations with the integration of cjj after the integration of is performed. 

Figure [5] shows the partial trace method for Eq. (j3-l|) . Figure a) shows the 
figure before the integrations are performed. Figure EJb) shows that the integration 
of the spin i is performed. Figure [2|^c) shows that the integrations of the spins i and 
j are performed. 



(a) 
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(b) 




Fig. 2. The partial trace method for Eq. (|3-l|l . (a) the figure before the integrations are performed, 
(b) the integration of the spin i is performed, (c) the integrations of the spins i and j are 
performed. 

By this partial trace method for Eq. (|3-ip . we are able to perform the integration 
of the partition fmiction Z. We obtain 

Wx,y} 

By using Eq. (j3-2p . we obtain 

By using Eq. (j3-3|) . we obtain the energy E as 

d 



E 



dp 



\n Z = -{L - if J tsx\hK . 



(3-2) 
(3-3) 

(3-4) 



The energy per site in the thermodynamic limit is obtained as — Jtanhi^. This 
result agrees with the previous resultP except for a constant and the coefficient of 
J. By using Eq. (|3-4p . we obtain the specific heat C as 



a/32 



C = kB/^'^TT^lnZ = ksK'^iL-lfsech'^K. 



(3-5) 
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By using Eq. ()3-3p . we obtain the free energy F as 

F = -^lnZ = -^[2(L-l) In2 + (L-1)2 ln(e^ + e"^)] . (3-6) 
p K 



By using Eqs. ()3-4p and ()3-6p . we obtain the entropy S as 

S=^iE-F) 

= 2{L - l)kB ln2 + A:b(L - 1)^ [ln(e^ + e"^) - KiauhK] . (3-7) 

We define the ground-state entropy per site in the thermodynamic limit as s'^^\ We 
obtain 

= lim hm — = . (3-8) 

K^ooN-^ooN 

The number of ground-state degeneracy, 17, is obtained as 

Q = hm e^/^B = 22(-^-i) . (3-9) 

K^oo 

Therefore, the number of ground-state degeneracy is not one, although s^^^ = 0. 

We calculate the thermal average of a spin, {(Jx,y)T- At first, we perform the 
integrations of all the spins except for the focused spin ax,y by applying the partial 
trace method for Eq. (j3-ip . Therefore, we obtain 



Ky>r = 0. (3-10) 

The present result shows {(Jx,y)T = even when T = 0. 

We calculate the two-point correlation function {(Jx,y(^x' ,y')T- At first, we perform 
the integrations of all the spins except for the focused spins ax,y and Ox'^y' by applying 
the partial trace method for Eq. ()3-ip . Therefore, we obtain 

{o'x,yO'x',y')T = ^x,x'^y,y' , (3-11) 

where 5xy is the Kronecker delta. 

When calculations of more than three-spin correlations, we perform the integra- 
tion of the spin at, which is one of the focused spins for thermal average, as 



f^i f^j Tfc 

(e-^ -e^)J2J2 <^j'^kA{a,,ak, {ai}) . (3-12) 



If A(aj,ak, {o'l}) is a Boltzmann factor for aj, and {cri}, we are also able to perform 
the integrations of aj and cjfc by using similar integrations with the integration of 
(Tj after the integration of is performed. The partial trace methods for Eqs. (j3-ip 
and (|3-12p are the partial trace methods that we propose in this article. 
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By applying the partial trace methods for Eqs. ()3-ip and ()3-12p . we obtain 

Wx,yCrx+i,y(^x,y+i)T = - tanh K . (3-13) 

Similarly, we obtain 

{(^x,yCrx+l,yC''x,y+l(^x',y'(^x'+l,y"^x',y'+l)T 

= (tanhi^)^ + [1 - (tanhi^)2](5^y5yy . (3-14) 

We calculate the three-point correlation function and the four-point correlation 
function. The three-point correlation function has been already derived by a different 
method)^ and the four-point correlation function has not been derived. We are able 
to consider the inside surrounded by the focused spins for thermal average and the 
outside of the spins separately. By using two functions vloutside and ^inside, we define 

^outside Ansido = ( CFi,j)TZ , (3-15) 

where {(i,j)} is the set of coordinates of the focused spins for thermal average, 
^outside is a Weight for outside of the focused spins, and ^inside is a weight for inside 
surrounded by the focused spins. We define the number of spins for inside surrounded 
by the focused spins as n, and define the number of plackets for inside surrounded 
by the focused spins as np. By applying the partial trace method for Eq. (j3-ip . we 
obtain 

^outside = 2^-^^-'^+"^(e^ + e-^)^^-'^^ . (3-16) 

We define the number of applications of the partial trace method of Eq. (j3-ip for 
inside surrounded by the focused spins as ni , and define the number of applications 
of the partial trace method of Eq. ()3-12p for inside surrounded by the focused spins 
as 722. There is a relation, written as 

np = ni + 71-2 . (3-17) 

If there is a correlation for the focused spins, we obtain ^inside as 

Ain,ide = 2"^"i-"2(g/^ + e'^y^ie-^ - e^T^ . (3-18) 

If there is no correlation for the focused spins, we obtain ^mside as 

Anside = 0. (3-19) 

At first, we calculate ^outside- The calculation of ^outside is started from one 
of the spins on th.G free bounda-riGS. -'^outside is calculated by applying the partial 
trace method for Eq. ()3-ip . Next, we calculate inside- The calculation of A\aside is 
started from one of the focused spins for thermal average, ^inside is calculated by 
applying the partial trace methods for Eqs. ()3-ip and p-12p . In a certain case, the 
partial trace method for Eq. (j3-ip is not applied to the calculation of j4inside5 while 
the partial trace method for Eq. (|3-12p is applied to the calculation of vlmside- We 
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find the spin correlations by examining the equations after the partial trace method 
for Eq. (j3-12p is applied. 

For the three-point correlation function, we generally obtain 

{o'x,y(^x',y'<^x",y")T = 0. (3-20) 

On the other hand, when {x',y') and {x",y") have a special relation, there is 

a correlation. When coordinates of three spins are {x, y), {x + 2^,y), (x, y + 2^) and 
k = 0,1, 2,..., there are correlations. Then, by applying the partial trace methods 
for Eqs. and we obtain n = {2^-^ + 1)(2'= + 1), m = 2^-^ {2^ + 1) - 3^ 

and n2 = 3^=. By using Eqs. ([321), (1345)1 - ([3481) . we obtain 

(o"x,j/0-^+2fc,yO-x,?;+2fc)T = -(tanhK)3\ A; = 0, 1, 2, . . . . (3-21) 

This result agrees with the previous result obtained in Ref . 3]) • 
For the four-point correlation function, we generally obtain 

{(^x,yO'x',y'(^x",y"0'x"',y"')T = 0. (3-22) 

On the other hand, when (x, y), (x', y'), (x", y") and (x'", y'") have a special relation, 
there is a correlation. When coordinates of four spins are (x,y), (x + 3 • 2^,y), 
(x,y + 3 • 2^), (x + 2'=,y + 2^) and k = 0,1,2,..., there are correlations. Then, 
by applying the partial trace methods for Eqs. (|3-ip and ()3T2p . we obtain n = 
(3 • 2^-^ + 1)(3 • 2^^ + 1), ni = 3 • 2'=-i(3 • 2^= + 1) - 6 • 3'' and na = 6 • 3''. By using 
Eqs. ([321), (I3T5H - (l34H|) . we obtain 

(cri^,i/0"z+3-2fc,j/Crx,y+3-2'=crx+2*,y+2fe)r = (tanhi^)*^'^ , A; = 0, 1, 2, . . . . (3-23) 

For the four-spin correlations, there are positive correlations. These positive corre- 
lations for spins have not been mentioned in the previous articles. 

From the solutions of the specific heat and the correlation functions, i.e., from 
Eqs. ()3-5p . ()3Tip . ()3-20p - ()3-23p . the transition temperature of this model is zero, 
although {(Jx,y)T = even when T = 0. This result, that the transition temperature 
is zero, is concluded from the thermal dependence of {(Jx,y)T in Ref. |^. Note that, 
in Ref. |4j), {ax,y)T = — 1 when T = and {(Jx.y)T = when T > are obtained. 

We investigate the spin configuration on a placket. We define the probabil- 
ity that the spin configuration {ax,y,crx+i^y,crx^y+i) on a placket is (— 1,— 1,— 1) as 
Px,yi-'i-, -1, -1)- -Px,y(-1, -1, -1) is given by 

Px,y{-^, -1, -1) = ^(1 - C7x,y){l - C7x+l,y){l " a,,y+i) . (3-24) 

We obtain 

2K 

(P.,,(-l,-l,-l))T = ^^^27r^- (3-25) 

We define the probability that the spin configuration {(Tx,y, ax+i,y, ax,y+i) on a placket 
is (-1, 1, 1) as Px,y{-1, 1, 1). Px,y(-1, 1, 1) is given by 

Px,y{-1, 1, 1) = -^(1 - o-x,s/)(l + cr^+i,y)(l + o-^.y+i) • (3-26) 
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We obtain 

2K 

(P.,,(-1,1,1))t = ^P^. (3-27) 

Similarly, we obtain 

2K 

{P^Jl, -1, = 1, -1))t = ^ . (3-28) 

When K — oo, {Px,y)T becomes the probability in the ground states. Therefore, 
in the ground states, the spin configurations (—1, —1, —1), (—1, 1, 1), (1, —1, 1) and 
(1, 1, —1) emerge with each l-in-4 chance. 

§4. Discussion 

There are disagreements between a part of the present results and the part of 
the previous results. We discuss the disagreements here. 

The present result shows that the number of ground-state degeneracy is 2'^^^~'^\ 
On the other hand, in Ref. [Tj), it has been pointed out that the number of ground- 
state degeneracy is one. 

The present result shows {(t)t = even when T = 0. In addition, the probability 
that the spin configuration {ax,y,(yx+i,y,(yx,y+i) on a placket is (— 1,— 1,— 1) is 1/4 
when r = 0. On the other hand, in Refs. it has been pointed out that 

{(t)t = — 1 when T = 0. 

This model has a special property. This model has been exactly solved in a 
special case that, under periodic boundary conditions, the lattice length is a power 
of two along at least one dimension, and has not been exactly solved on lattices of 
any length, so that there is a possibility that a part of the results on lattices of any 
length is different from the part of the previous results. From the results on lattices 
of any length under free boundary conditions, we see that there are the differences. 
Especially, the results for the ground states are different from the previous results 
as described above. 

These disagreements between the present results and the previous results would 
not seriously affect the dynamics. The dynamics are decided by the energy changes 
of spin-flips,'^ and the energy changes of spin-flips are not mentioned in this article, 
but the energy changes of spin-flips are the same with the previous results. The 
number of ground-state degeneracy is different from that of the previous result, but 
the ground state entropy in the thermodynamic limit is the same with the previous 
result, so that the disagreement for ground-state degeneracy would not seriously 
affect the dynamics. Since the dynamics are decided by the energy changes of spin- 
flips, the effect from the different boundary conditions would not seriously affect the 
dynamics if the system size is large enough. 

§5. Summary 

A spin model, which exhibits glassy dynamics and has no randomness for ex- 
change interactions, has been proposed by Newman and Moore. We proposed partial 
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trace methods for obtaining static exact solutions of this model under free boundary 
conditions, and showed some static exact solutions. The partial trace methods for 
Eqs. and (j3-12|) are the partial trace methods that we proposed in this article. 
We obtained that, from the exact solutions of the specific heat and the correlation 
functions, the transition temperature is zero, although the thermal average of a spin 
is zero when the temperature is zero. We investigated the ground states. In the 
ground states, a part of the present results especially disagreed with the part of the 
previous results. We discussed the disagreements, and concluded that the disagree- 
ments would not seriously affect the dynamics. 
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